Visualizing elusive phase transitions with geometric entanglement 
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We show that by examining the global geometric entanglement it is possible to identify "elu- 
sive" or hard to detect quantum phase transitions. We analyze several one-dimensional quantum 
spin chains and demonstrate the existence of non-analyticities in the geometric entanglement, in 
particular across a Kosterlitz-Thouless transition and across a transition for a gapped deformed 
Ameck-Kennedy-Lieb-Tasaki chain. The observed non-analyticities can be understood and classi- 
fied in connection to the nature of the transitions, and are in sharp contrast to the analytic behavior 
of all the two-body reduced density operators and their derived entanglement measures. 



I. INTRODUCTION 

The effect of interactions in many-body systems gives 
rise to striking collective phenomena 1 . Phase transitions, 
both classical and quantum, are the archetypical exam- 
ple. Across such a transition, collective properties of the 
system undergo abrupt changes that can sometimes be 
related to non-analytic behavior of the free energy. This 
observation was at the basis of the first historical attempt 
to classify phase transitions by Ehrenfest, according to 
the order of the non-analyticity involved. Modern classi- 
fication schemes have refined this idea in order to include 
new types of transitions 1,2 . 

Across quantum phase transitions (QPT), one expects 
that the ground-state wavefunction undergoes drastic 
changes and hence manifests this behavior via physi- 
cal quantities such as correlations. Recently there has 
been a significant effort towards exploring the relation 
between the revived quantum-mechanical entanglement 
and QPT to complement traditional approaches. For in- 
stance, important scaling properties have been found for 
the entanglement entropy and single-copy entanglement 
between macroscopic regions in various contexts J , includ- 
ing the connection to the central charge. A different ap- 
proach has been the use of entanglement between individ- 
ual constituents, such as the two-qubit concurrence and 
other correlation-based measures . In particular, con- 
currence was demonstrated to display singularity across 
QPT . It was later recognized that such non-analytic 
behavior originates in the two-body reduced density ma- 
trices and is linked to the non-analyticity in the ground- 
state energy (the so-called "generalized Hohenberg-Kohn 
Theorem" ) ' . Also, a similar approach (but not origi- 
nated from entanglement) called fidelity measure, which 
employs the overlap between two ground states, has been 
successful in identifying QPT 11 . 

According to the above picture, it is possible to de- 
tect finite-order transitions just by examining the non- 
analyticities of two-body entanglement measures. How- 
ever, one encounters difficulty with other types of tran- 
sitions. For instance, in oo-order transitions, such as 
Kosterlitz-Thouless (KT), the ground-state energy and 
its derivatives are analytic, as well as all correlation func- 
tions, such as two-body observables. This is the case of, 



e.g., the spin-1/2 XXZ chain near the antiferromagnetic 
Heisenberg point. A further example, not of the KT type, 
is a transition occurring in a deformed Afncck-Kennedy- 
Lieb-Tasaki (AKLT) chain introduced by Verstraete et 
al. in Ref. N , where the existence of a diverging entan- 
glement length scale in the system remains undetected 
by any correlation functions of the ground state, as the 
system is always gapped and the correlation length re- 
mains finite. The complex nature of these transitions 
makes them elusive and undetectable by all the above 
entanglement approaches (as well as the fidelity suscep- 
tibility measure) " 1 , and previous investigations indicate 
that they may be better understood in terms of global 
quantities 1 1-13 . 

Here we provide a new perspective along this direction, 
and show that for ID quantum many-body systems the 
global geometric entanglement can be used to successfully 
detect QPT, including finite-order and the above elusive 
ones. The geometric entanglement (GE) 1 ' has previously 
been shown to exhibit divergence consistent with a scal- 
ing hypothesis 15 ' 18 , and has also been related to the cen- 
tral charge of the underlying conformal theories at criti- 
cality 11 . Moreover, its finite-size corrections at criticality 
are also governed by conformal symmetry . In this con- 
text, the aim of the present work is to show that even 
when all correlation functions remain analytic, the GE 
is still able to display singularity across transitions. We 
shall also elaborate on the connection between the ori- 
gin of these singularities and the nature of the observed 
transitions. 

The structure of this paper is as follows: in Sec. II we 
review briefly the basics on the global geometric entangle- 
ment. In Sec. Ill we show our results for a variety of ID 
systems, namely the spin-1/2 Ising model in transverse 
and longitudinal fields, the spin-1/2 XXZ model, and the 
deformed AKLT model. Sec. IV offers a discussion of 
the results focusing on two aspects: first, the connection 
between the observed singularities for the GE and the 
nature of the phase transitions, and second, a compari- 
son of the performance to detect QPTs between the GE 
and other entanglement-related quantities. Finally, Sec. 
V contains the conclusions. 
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II. GLOBAL GEOMETRIC ENTANGLEMENT 
PER SITE 

In this section we briefly remind the basics of the global 
geometric entanglement. To characterize global entan- 
glement, consider a general, TV-partite, normalized pure 
state \V) £H = <S)f =1 H [i] , where H H is the Hilbert 
space of party i. For a spin system each party could be 
a single spin, but could also be a block of contiguous 
spins ' . Our scheme involves considering how well an 
entangled state can be approximated by some unentan- 
gled (normalized) state: |3>) = ®£Li \4$), motivated by 
mean-field theory. For quantum spin systems, the mean- 
field scheme attempts to find the best product state |$) 
minimizing the Hamiltonian H. Here, we aim to find the 
best mean-field approximation to the ground state 
The proximity of \^f) to |<&) is captured by their over- 
lap; the entanglement of is revealed by the maximal 
overlap : 

A max (*) = max|($|*)|. (1) 

A is thus related to the best mean-field energy of the 
"reduced" Hamiltonian 

H Ied = -\V}{H>\, (2) 

as the closest product state |3>*) is the best mean-field 
state such that 

min<*|ft rcd |*) = ($*|ff red |$*) = -A max (*) 2 . (3) 

It makes sense to quantify the entanglement via the fol- 
lowing extensive quantity 15,18 (analogous to the relation 
between the free energy and the partition function): 

£(#) = - log 2 AL^*), (4) 

where we have taken the base-2 logarithm. GE gives zero 
for unentangled states and is a measure of how difficult it 
is to approximate a given state (in particular the ground 
state) by mean-field states, or equivalently a measure of 
unfactorizability. To deal with large systems we define 
the thermodynamic entanglement density £ and its finite- 
size version £ jy by 

£ = lim £ N , £ N = ^1. (5) 

This is the quantity that will be of interest in this paper. 

III. VISUALIZING DIFFERENT TYPES OF 
TRANSITIONS 

In this section we provide an analysis of different 
ID quantum spin systems undergoing different types of 
QPTs, from the point of view of the GE. Specifically, 
we focus on the spin-1/2 Ising model in transverse and 
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FIG. 1: (color online) £ for the Ising model obtained with 
MPS vs. the transverse field h for several values of the longi- 
tudinal field A (left) and vs. the longitudinal field A for several 
values of the transverse field h (right). The insets show the 
derivatives with respect to h and A. The derivative d£/dh in 
the left panel corresponds to the line for A = 0. 

longitudinal fields, the spin-1/2 XXZ model, and the de- 
formed AKLT model. Let us mention that the global ge- 
ometric entanglement per site £ has already been applied 
to ground states of ID models across different types of 
phase transitions 14,15 ' 18 . Our analysis here complements 
those from previous studies by offering new results in 
more exotic situations. 

A. Spin-1/2 Ising model 

For comparative purposes, we first revisit the spin- 
1/2 quantum Ising model in transverse and longitudinal 
fields, 

H=-J2(aMo-^+ha®+\a®), (6) 

i 

where h (A) is the transverse (longitudinal) field, and 
a a is the a-th Pauli matrix at site k. This Hamilto- 
nian has a Z2 symmetry-breaking second-order quantum 
phase transition at h* = ±1 and A = 0, whereas at fixed 
\h\ < 1 it has a first-order discontinuous transition at 
A* = due to a crossing of energy levels. Namely, the 
phase diagram is a first-order line terminated by second- 
order points. For this model, we employ the infinite- 
TEBD algorithm 1 '' to find a MPS approximation to the 
ground state in the thermodynamic limit. Then, GE is 
readily obtained from the MPS state by maximizing the 
overlap (1) with standard optimization. 

Results for the ground state of the quantum Ising 
model in Eq. (6) are shown in Fig. 1. On the left panel we 
extend the results from Ref. )s for GE across the quantum 
phase transition as a function of the transverse field h for 
different values of A. We have checked that our MPS re- 
sults for A = reproduce accurately the exact solution 
in Ref. (less than 1% of relative error). Notice that, at 
A = 0, £ is smooth across the second order phase transi- 
tion with a peak slightly after the quantum critical point 
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FIG. 2: (color online) £n for the XXZ model in zero field 
vs. the anisotropy parameter A for system sizes N = 
10, 12, 14, 16, 18, 20, 24 and 26. The dashed lines correspond 
to the thermodynamic limit £, obtained by fitting the finite- 
size scaling law in £m{A) ~ £(A) + b(A)/N (upper dashed 
line) and by the infinite MPS method (lower line). We also 
indicate the closest product state |<3>) on each side. 



FIG. 3: (color online) £ N for the spin-1/2 XXZ model 
in zero field, as a function of 1/N for system sizes N = 
10, 12, 14, 16, 18, 20, 24 and 26. The data correspond to A = 
0.5, 1 and 1.5. The dashed lines are our best fits to Eq. (10). 
The values extrapolated to the thermodynamic limit 1 /N — > 
correspond to those of the dashed line plotted in Fig. (2). 



B. Spin-1/2 XXZ model 



(around h ~ 1.13). The derivative, however, is divergent 
at the quantum critical point ft = h* = 1, as shown in 
the inset, and obeys the critical scaling law 



-(A = 0,ft)~--log 2 |ft-l| (7) 

for | ft — 1| <C 1 ■ Also, as can be easily inferred from 
Fig. 1, our MPS results prove that this derivative is 
smooth for A ^ 0, as there is no transition. The be- 
havior of £ is rather different across the line of the first- 
order transition as a function of the longitudinal field A, 
for which we give our results in the right panel of Fig. 1. 
There, we see that £ has a kink (thus being non-analytic) 
as a function of A at the first-order (discontinuous) phase 
transition point A = A* = for ft ^ 0, 1. At the second 
order phase transition point A = 0, ft = 1 our MPS re- 
sults are compatible with a logarithmic divergence of the 
derivative 



— (A,ft = l)~-alog 2 |A| (8) 



for |A| <C 1, with a ~ —7.5(1). Notice that £ is symmet- 
ric around A = A* since at this point the Hamiltonian is 
self-dual under the duality transformation A — > —A. The 
observed non-analiticity at A = and /i ^ 0,1 can be 
understood as a consequence of a global change in the 
ground-state wavefunction, characteristic of first-order 
transitions with a crossing of ground state energy levels. 
What is more intriguing is that similar non-analytical be- 
haviors in £ are also found in other types of transitions, 
even continuous ones, as we will see in what follows. 



We now consider a system with an elusive phase tran- 
sition, i.e., the ID spin-1/2 XXZ model, 

H = J2 + + A ^ ] 4 i+11 + haf) , 

i 

(9) 

where A is an anisotropy parameter and h a magnetic 
field. 

Let us first study the case of zero field (ft = 0). In 
this regime, this model is critical for A € (—1,1], with 
a KT quantum phase transition at the Heisenberg point 
A* = l 21 . Within this setting, first we do an exact di- 
agonalization of H for sizes up to 26 spins and find the 
geometric entanglement, followed by a finite-size scaling 
and extrapolation to the thermodynamic limit. Then 
we compare this value with that obtained by using the 
MPS method for infinite systems, as used for the quan- 
tum Ising model. In turn, this allows us to further vali- 
date the consistency of our numerical methods. 

In Fig. 2 we show the results for the spin-1/2 XXZ 
model (9) in zero field. We see that the global geometric 
entanglement per site £n for finite size N already displays 
a pronounced kink at the KT quantum critical point A = 
A* = 1. As observed in the figure, £n seems to converge 
fast towards a thermodynamic value as N increases. Our 
finite-size scaling analysis indicates a scaling law 

£ N (A)~£(A)+ b -^, (10) 

in good agreement with the ones proposed in Ref. 15 , see 
Fig. 3. 

We have done this scaling analysis for all the computed 
values of A and obtained an estimation of the thermo- 
dynamic quantity £, shown in Fig. 2, together with the 
infinite MPS estimation and the finite-size data. The val- 
ues of £ estimated by both methods agree within 1% of 
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FIG. 4: Phase diagram of the XXZ model from Eq. (9) in the 
(h, A) plane. The dots correspond to our estimation observing 
the discontinuities of the GE using MPS methods for infinite 
systems, and the lines are the exact phase boundaries. Phases 
A and C are gapped, whereas phase B is gapless. Phases 
A and B are separated by a line of first order transitions, 
whereas phases B and C are separated by a line of second 
order transitions at ft > that ends in a Kosterlitz-Thouless 
(KT) transition at h = 0. 



relative error, which validates the consistency of our dif- 
ferent approaches. The kink in £ jy at A = 1 is obviously 
present in the thermodynamic limit N — > oo. This is a 
remarkable result, as for this KT transition all the two- 
body observables and all their derivatives are analytic, 
and this means that entanglement measures that only de- 
pend on two-body reduced density operators, such as the 
concurrence and the spin-spin negativity, will not exhibit 
any singularity at all. Our results also indicate that this 
kink is due to a sudden change in the product state that 
maximizes the overlap in Eq. (1): for A < 1 the closest 

product state |$) is | • • H !-—•••) (as well as those from 

rotating | • • • -\ 1 — ...) around z axis, due to SO(2) 

symmetry), whereas for A > 1 it is | . . . 0101 . . . ) (where 
|±) = (|0) ± |l))/\/2, and |0) and |1) are the eigenstates 
of <7 Z ). At the isotropic point A = 1, either of the two 
product states is equally good, due to the SU(2) symme- 
try. The observed kink in the global geometric entangle- 
ment evidences the existence of the KT transition, and 
its similitude with the non-analytical behavior found in 
discontinuous phase transitions (see Fig. 1) supports the 
fact that there is a sudden and global change in the struc- 
ture of the ground-state wavefunction . Notice, though, 
that according to the standard classification, the phase 
transition in this model is continuous. 

Let us now consider the case of nonzero field (h ^ 0) in 
the Hamiltonian. In this generic case, by looking at the 
discontinuities in the GE, we find out that it is also pos- 
sible to get both a qualitative and quantitative picture 
of the correct phase diagram for this model in the (h, A) 
plane. Our results for the phase diagram estimated using 
the GE obtained from MPS methods for infinite systems 
are shown in Fig. 4. As seen in the plot, there is a good 
quantitative agreement between the phase boundaries es- 
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FIG. 5: (color online) £ for blocks of 2 spins for the deformed 
AKLT model vs. the deformation parameter /j,. The right 
inset shows a close-up around the AKLT point fx* = 0, and 
the left inset shows the coefficient a(n) associated with the 
closet product state |$(^i)) (see text). These results are exact. 

timated with the GE and the exact ones (computed by 
Bethe ansatz). Thus, we see that the GE is able to re- 
produce within good accuracy the correct properties of 
the phase diagram of the model. 

C. Deformed AKLT model 

Finally, we consider the deformed AKLT modeP, 
H = ]TAf +1 ] (11) 

i 

x [i,i+i] = ((sWj-^sMl^f^])-^^]), 

where S^ 1 = iW + sinh (p)S^ ] + (cosh (/x) - l)(sf ] ) 2 , and 

*ak£t = ■ + Us® ■ S^) 2 + \ (12) 

is the usual AKLT two-body term 1 ' 1 , with S$ the a- 
th component of the spin-1 operator S^. The ground 
state undergoes a transition at the AKLT point \x* = 0, 
with diverging entanglement length and finite correlation 
length. For this model, we use the exact MPS represen- 
tation of the ground state from Ref. and then extract 
from it the geometric entanglement per site in the ther- 
modynamic limit by using the exact MPS techniques de- 
veloped in the second paper of Ref. 1 '. For convenience of 
the calculation, now we choose each party in Eq. (5) to 
be composed by a block of two contiguous spins, so that 
£ now refers to the geometric entanglement per block of 
two spins in the thermodynamic limit . 

From Fig. 5 we see that GE £ has a pronounced kink at 
the AKLT point fi = fi* — in the thermodynamic limit, 
similar to that in the KT transition of the XXZ model 
and the first-order transition of the Ising spin chain in a 
longitudinal field. This similitude supports again the idea 
of a sudden global change in the ground-state wavefunc- 
tion. However, notice that as explained in Ref. v , here the 
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system is always gapped and the correlation length of the 
ground state of the system is always smooth and remains 
finite for this transition. Thus, this sort of transition 
does not even exist according to the standard criteria, 
and two-body correlation functions are unable to detect 
the observed non-analyticity. However, the entanglement 
length diverges at the AKLT point . Remarkably, we see 
here that GE is also successful in identifying the exis- 
tence of this transition in the ground state of the system. 
We also determine the closest product state 



l*(A0>= C(/x)U(/i)|0,0) + i|x( M )) 



(13) 



with 



|z(/i<0)> = |-1,-1>, \x(jm>0)) = \1,1), (14) 

C(/i) a normalization constant, a(/x) a real positive coef- 
ficient (see Fig. 5), and | — 1), |0) and |1) the eigenstates 
of the spin-1 operator S z . 



IV. DISCUSSION 

The results obtained in the previous section prove the 
usefulness of the GE to detect phase transitions of many 
different kinds, including those that seem difficult to de- 
tect with alternative methods. In this section we dis- 
cuss a number of important considerations that can be 
observed from our results, namely (i) how do the singu- 
larities in GE connect to the nature of the transitions, 
and (ii) how does the GE compare to other alternative 
measures in efficiency of calculations and visualization of 
results. 



state | , and (ii) parameters r * M , which denote the vec- 
tors that maximize the overlap. 

In all the examples that we have examined, we can 
classify the origin of the singularity due to (i) or (ii) or 
both. Let us summarize: 



1. For the transverse Ising model, which has a stan- 
dard second-order quantum critical point, (i) cor- 
relation functions C's are singular but (ii) optimal 
parameters r*'s are not singular. This explains 
the similar behavior between the GE and the so- 
called concurrence measure of entanglement, which 
depends on correlation functions. 

2. For the longitudinal Ising model, which has a stan- 
dard first-order transition, both (i) and (ii) are sin- 
gular, as the transition is first-order. 

3. For the XXZ model at zero field, the transition is 
oo-order, therefore (i) correlations C's are not sin- 
gular, but (ii) the parameters r*'s of the optimal 
local states are singular. It is this second point 
the one that helps to signify certain non-analytic 
change in the wavefunction and thus identifies the 
transition. 

4. For the deformed AKLT model, (i) correlations C's 
are not singular, since the correlation length is fi- 
nite, but (ii) r*'s are singular. Similar to XXZ, 
it is this second point the one that detects non- 
analyticity in the wavefunction across the transi- 
tion. 



A. Singularities in GE and the nature of the 
transitions 

To illustrate the nature of the observed singularities 
in the GE, consider for simplicity spin-1/2 systems and 
rewrite the A-spin product state |4>) via 

|$)<$| = | I(lM+f*] .5*1), (15) 

i—l 2 

where the unit vectors fM 's represent directions of local 
spins. The overlap K^l^)! 2 can be expressed (by ex- 
panding the above product) as a linear combination of 
all correlations w.r.t. 

2 N m n 2 = l + E^-(* ll] )+ E ^rf{a^f) 

+ E r®rfr!?H*®<rfaW) + --.. (16) 

This can be easily generalized to systems of higher spins. 
Therefore, it is seen that a singularity of the entangle- 
ment can come from two types of sources: (i) correla- 
tion functions, C^g*'"'} = (ca tr^V^ . . . ) for the ground 



B. Comparison to correlation functions 

Let us now discuss the relevance of the GE as com- 
pared to other approaches based on correlation functions 
to study phase transitions. One could be tempted to af- 
firm that any phase transition, if it exists, can in princi- 
ple be detected by measuring all the possible correlation 
functions of the system, and that therefore the GE of- 
fers no true extra information and is not useful to study 
phase transitions. 

We argue here that this approach may not actually 
apply. To see this, first notice that it is impractical to 
exhaust all possible correlations to find out if there is any 
singularity in a given system. And moreover, there exist 
example Hamiltonians where all the ground-state corre- 
lations are well-behaved and no singularity can be found 
(e.g. XXZ model at h = 0). Entanglement measures that 
depend on correlation functions, such as the concurrence, 
also inherit this analytical behavior. Therefore, a quan- 
tity which includes all the possible correlation functions 
in a single quantity is potentially very useful as one then 
would need to examine this single quantity to see if any 
singularity exists in the correlations. As shown above, 
one has that (i) the GE is actually such a quantity since 
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it can be expressed in terms of a combination of all pos- 
sible correlation functions (general k-point correlations), 
and (ii) there are additional quantities (e.g. the vectors 
that characterize the best local product state) that also 
assist the examination of singularities. In all our exam- 
ples examined in the paper, we can classify the origin of 
the singularity is due to (i) or (ii) or both. In short, GE 
is a simple and meaningful quantity that provides infor- 
mation that cannot be codified in any (local) correlation 
function, and the approach is clearly more efficient than 
calculating all possible correlation functions and exam- 
ining them one by one. 

C. Comparison to other measures: localizable 
entanglement, entropy and fidelity 

Let us now briefly discuss how the GE compares to 
other entanglement-related quantities in detecting phase 
transitions. We focus on the localizable entanglement, 
the entanglement entropy, and the ground-state fidelity. 
Notice that of all these quantities, the fidelity is not a 
measure of entanglement by itself. We remark, however, 
that this does not diminish its usefulness in studying 
phase transitions. 

Let us first considered the appearance of singular be- 
haviors across QPTs. Quite importantly, for the two 
elusive models studied in this work (XXZ and deformed 
AKLT), measures such as the entanglement entropy and 
fidelity measures 1 (which we also analyzed for deformed 
AKLT - results not shown -) fail to show any singularity 
across the transitions. Furthermore, if one considers a 
derived quantity from fidelity, called fidelity susceptibil- 
ity, it can be shown that for KT transition it does not 
show any singularity . The localizable entanglement, 
however, shows a singular behavior in these two transi- 
tions as well 8,12 , in a way similar to the one observed 
with the GE. Notice, though, that the GE may be easier 
to compute than the localizable entanglement in many 
cases, as we argue below. 

Next, let us consider the efficiency in the calculation 
of the different measures. For certain exactly solvable 
models, the GE can be calculated essentially analytically. 
Furthermore, for non-solvable models, with existing nu- 
merical techniques based on tensor networks such as MPS 
or PEPS Ju it is rather straightforward to compute the 
GE. 

In fact, the geometric entanglement is probably one 
of the multipartite entanglement measures that is easier 
to calculate, while most of the other known multipartite 
entanglement measures end up being rather hard to com- 
pute. For example, in the framework of MPS, computing 
GE is not harder than computing the ground state. Once 
we have a MPS approximation of the ground state, it is 



numerically easy to calculate the GE, i.e., we take bond 
dimension 1 in an MPS that minimizes the energy of 
the Hamiltonian W = -|*)(*|, where |f) is the MPS 
ground state approximation of the original Hamiltonian 
under study. Comparing to other approaches, we note 
that the calculation of the fidelity susceptibility also re- 
quires the knowledge of ground states, but is perhaps 
more inefficient than our approach from a computational 
point of view, as it would require the overlap between two 
MPS instead of an MPS and a product state. Regarding 
the entanglement entropy, its calculation requires com- 
puting the reduced density matrix of a block of finite 
size together with its spectrum, which can not always be 
done efficiently (especially for systems in more than one 
dimension). Finally, regarding the localizable entangle- 
ment, there is in general the necessity to maximize over 
all possible local measurements (not necessary projective 
measurements), which makes it the most difficult calcu- 
lation of all the ones discussed so far for a generic model. 



V. CONCLUSIONS 

Here we have shown that the global geometric entan- 
glement can be used to successfully detect conventional 
and elusive phase transitions, such as the ones for the ID 
spin-f/2 Ising, XXZ and deformed spin-1 AKLT models. 
We have also clarified the connection between the nature 
of the observed transitions and the fact that the geo- 
metric entanglement exhibits singularities whereas other 
entanglement measures do not. Thus, we have demon- 
strated that the GE can be used to detect elusive tran- 
sitions for which other conventional methods and other 
entanglement measures (including the fidelity suscepti- 
bility between ground states) fail. 

All in all, we believe that the GE may provide com- 
plementary information about the complicated ground 
states of quantum many-body systems to the one offered 
by alternative measures such as e.g. correlation func- 
tions, entanglement entropy, localizable entanglement 
and ground state fidelity. However, there is still lack of 
extensive study of the relations amongst all these meth- 
ods. Further study in this direction would help to clar- 
ify the complex nature of the ground state of strongly- 
correlated systems. 
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